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O . 1 INTRODUCTION 
t-H . 

In the restricted three-body problem there are different domains of 
stable motion associated to co-orbital motion. Each can be clas- 
, sified according to the center of libration of the critical argument, 
$_( ' a = A — A', where A denotes the mean longitude of the minor body 
and A' the same variable for the disturbing planet. These types of 
motion are known as: (i) Tadpole orbits, corresponding to a libra- 
tion of a around L4 or L5; (ii) Horseshoe orbits, where motion 
occurs around a = 180° and encompasses both equilateral La- 
grangian Points, and (Hi) quasi-satellite (QS) orbits, where a oscil- 
lates around zero. 

The term "quasi-satellites" was originally introduced by 
Mikkola and Innanen (1997) and can be viewed as an extension 
of retrograde periodic orbits in the circular restricted three-body 
problem (e.g. Jackson 1913, Henon 1969). Although not present 
for circular orbits, they exist for moderate to high eccentricities of 
the particle. In a reference frame rotating with the planet, QS orbits 
circle the planet like a retrograde satellite, although at distances so 
large that the particle is not gravitationally bounded to the planetary 
mass (Mikkola et al 2006). 

The first object confirmed in a QS configuration was the 
asteroid 20021/S68 (Mikkola et al., 2004) with Venus as the 
host planet. The Earth has one temporary co-orbital object, 
(3753 Cruithne, Namouni et al. 1999), and one alternating 
horseshoe-QS object (2002AA29, Connors et al. 2002). The co- 
orbital asteroidal population in the inner Solar System was stud- 



ABSTRACT 

We study the stability regions and families of periodic orbits of two planets locked in a co- 
orbital configuration. We consider different ratios of planetary masses and orbital eccentric- 
ities, also we assume that both planets share the same orbital plane. Initially we perform 
numerical simulations over a grid of osculating initial conditions to map the regions of sta- 
ble/chaotic motion and identify equilibrium solutions. These results are later analyzed in more 
detail using a semi-analytical model. 

Apart from the well known quasi-satellite (QS) orbits and the classical equilibrium La- 
grangian points L4 and L5, we also find a new regime of asymmetric periodic solutions. For 
low eccentricities these are located at (AA, Aw) — (±60°, =f120°), where AA is the differ- 
ence in mean longitudes and Am is the difference in longitudes of pericenter. The position 
of these Anti-Lag rangian solutions changes with the mass ratio and the orbital eccentricities, 
and are found for eccentricities as high as ~ 0.7. 

Finally, we also applied a slow mass variation to one of the planets, and analyzed its effect 
on an initially asymmetric periodic orbit. We found that the resonant solution is preserved as 
long as the mass variation is adiabatic, with practically no change in the equilibrium values of 
the angles. 

Key words: celestial mechanics, planets and satellites: general. 



ied in Brasser et al. (2004) by numerical integrations. All QS or- 
bits appear to be temporary, escaping in timescales of the order of 
10 2 - 10 4 years 

Wiegert et al. (2000) numerically investigated the stability of 
QS orbits around the giant planets of the Solar System. Although 
no stable solutions were found for Jupiter and Saturn, some initial 
conditions around Uranus and Neptune lead to QS orbits that sur- 
vive for timescales of the order of 10 yr. It thus appears that a 
primordial population of such objects may still exist in the Solar 
System. Kortenkamp (2005) used N-body simulations to model the 
combined effects of solar nebula gas drag and gravitational scatter- 
ing of planetesimals by a protoplanet. He showed that a significant 
fraction of scattered planetesimals could become trapped into QS 
trajectories. It then seems plausible that this trapped-to-captured 
transition may be important not only for the origin of captured 
satellites but also for continued growth of protoplanets. 

At variance with these results, in the case of the general (non- 
restricted) three-body problem, although equilateral solutions and 
horseshoe orbits are well known, quasi-satellite configurations have 
only been studied very recently. Hadjidemetriou et al. (2009) per- 
formed a detailed study of periodic orbits in the 1/1 MMR for fic- 
titious planetary systems with different mass ratios. They found 
that stable QS solutions occur for a — AA = A2 — Ai = and 
Avj = VJ2— vj\ = 180° , where the subscripts identify each planet. 
Unstable trajectories were found at a = 180°, Avj — 0. Although 
at present there are no confirmed cases of exoplanets in quasi- 
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Figure 1. Results of numerical integrations of initial conditions in a grid in 
the (a, Aw) plane. Planetary masses were taken equal to mi = m,2 = 
rrt j U p, and initial semimajor axes equal to a± = a 2 = 1 AU. Regions of 
regular motion are shown in white, while the dashed regions correspond to 
chaotic and unstable trajectories. 



satellite configurations, Gozdziewski & Konacki (2006) found that 
the radial velocity curves of the HD82943 and HD 1283 11 planets 
could correspond to co-orbital motion in highly inclined orbits. Nu- 
merical simulations of both systems show QS trajectories, instead 
of Trojan orbits as initially believed. 

In the present work we aim to revisit the 1/1 mean-motion res- 
onance (MMR) in the planar planetary three-body problem, trying 
to identify possible domains of stable solutions and their location 
in the phase space. Section[2]presents several dynamical maps con- 
structed from numerical simulations for different initial conditions. 
These maps allow us to identify stable fixed points and periodic or- 
bits, as well as the domains of regular motions. In Section[3]we de- 
velop a semi-analytical model for co-orbital planets, which is then 
applied in Section [4] to calculate the families of stable periodic or- 
bits. In the same section we also present a brief study of the effects 
of an adiabatically slow mass variation in one of the planetary bod- 
ies. Finally, conclusions close the paper in Section|5] 



2 DYNAMICAL MAPS WITH EQUAL MASS PLANETS 

Consider two planets with masses mi and m2 in coplanar orbits 
around a star with mass mo = Mq . We will begin considering the 
case m,2 = mi, also other mass ratios will be discussed in later 
sections. Let a% denote the semimajor axes, e,: the eccentricities, 
Xi the mean longitudes and Wi the longitudes of pericenter. All 
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Figure 2. Semi-amplitude maps. The Left (Right) column shows the ampli- 
tude variation for a (Aro) in gray scale. Light domains correspond to near 
zero amplitude families, darker regions indicate oscillation amplitudes up 
to 90°, and dashed regions correspond to unstable orbits. Initial values of 
eccentricities are indicated in each panel. Color scale is indicated at bottom 
and ACR solutions are marked on the right panels. 



orbital elements considered in this paper are assumed astrocentric 
and osculating. Throughout this work, mi will be our "reference" 
planet: its mass will be fixed at one Jovian mass (mi = Mj up ) 
and the system scaled to initial condition a\ = 1 AU. The angular 
variables for co-orbital motion will then be defined as a = A2 — Ai 
and Auj — VJ2 — 

As pointed out by Hadjidemetriou et al. (2009), for equal mass 
planets the periodic orbits are such that are located at a\ = 02 and 
ei = 62. Accordingly, we fixed the semimajor axes and eccentrici- 
ties, and constructed a 100 x 100 grid of initial conditions varying 
both g and Au? between zero and 360°. Each point in the grid 
was then numerically integrated over 3000 orbital periods using a 
Bulirsch-Stoer based N-body code, and we calculated the averaged 
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MEGNO chaos indicator (Y) (Cincotta & Simo, 2000) to identify 
regions of regular or chaotic motion. Results are shown in Figure 
Q]for six values of the initial eccentricities e;; dashed regions cor- 
respond to unstable orbits while white was used to identify stable 
solutions. An analysis of these plots show the following character- 
istics: 

• For low initial eccentricities (e^ = 0.05) the maps show two 
disconnected strips of regular motion, corresponding to motion 
around a = ±60° and any value of Azu. 

• For moderate low to intermediate initial eccentricities (e^ = 
0.15 and e, = 0.30) the vertical strips of regular motion become 
thinner and slightly distorted. A new stable domain is now present, 
associated to QS orbits, and located around a = 0. 

• For high initial eccentricities (e; > 0.40) the domain of QS 
orbits increases and covers a significant portion of the plane of ini- 
tial conditions. Conversely, the distorted vertical strips shrink and 
each seems to break into two islands of stable motion. The smaller 
islands encompass equilateral solutions, although they almost dis- 
appear for ej = 0.70. The larger islands correspond to a different 
type of asymmetric solution, and their locations tend towards the 
center of the plots as the eccentricities increase. 

• Due to symmetry present in the dynamical system, the re- 
sults are invariant to transformations of the type (a, Azu) — > 
(—a, —Aw). In fact, since mi = rri2, both equilateral solutions 
are actually the same solution, since we can pass from one to the 
other just by redefining the reference planet. However, since later 
sections will discuss the case mi 7^ mi, we prefer to treat both 
equilateral solutions separately. 

Although MEGNO is a very efficient tool to identify chaotic 
motion, it is not suited to distinguish between different types of 
regular orbits (e.g. fixed points, periodic orbits, etc.). Sometimes 
this task is performed with a Fourier transform of the numerical 
data (e.g. Michtchenko et al. 2008ab); however, here we have cho- 
sen a different route. Starting from the output of each numerical 
simulation, we calculated the amplitudes of oscillation in each an- 
gular variable. Initial conditions with zero amplitude in a corre- 
spond to o-family periodic orbits of the co-orbital system, while 
solutions with zero amplitude in Azu will correspond to periodic 
orbits of the so-called Aw-family (see Michtchenko et al. 2008ab). 
Finally, stationary solutions of the averaged problem, identified as 
intersections of both families, may be thought as analogous to the 
apsidal corotation resonances (ACR) found in other mean-motion 
resonances (e.g. Beauge et al. 2003). The equilateral Lagrangian 
solutions will appear as ACR in these plots. 

The gray scale graphs in Figure [2] show values of the ampli- 
tudes in a (left) and Azu (right) for four of the plots shown in figure 
[TJ White regions represent initial conditions with semi-amplitudes 
smaller than 2°, as thus indicate the families of periodic orbits in 
each angle. Darker regions correspond to increasing amplitudes 
(up to 45°) and denote initial conditions with quasi-periodic mo- 
tion. The dashed areas are unstable solutions. Finally, it is worth- 
while mentioning that symmetric configurations may either corre- 
spond to an alignment (Azu = 0°) or an antialignment of the apses 
(Azu = ±180°) while asymmetric configurations have stationary 
values of Azu different from the above. 

For low eccentricities (e^ = 0.05) we observe four asymmet- 
ric ACR solutions. Two are the well known Lagrangian equilateral 
solutions located at (a, Azu) = (±60°, ±60°). By analogy with 
the restricted problem, we will denote them L4 and L5. As far as we 
know, the remaining two ACR have not been previously reported, 
and are located at approximately (a, Azu) = (±60°, =Fl20°). We 
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Table 1. Approximate location for the stable ACR solutions in the (cr, Aro) 
plane, calculated from the dynamical maps with ei = = 0.15. Forequal 
mass planets, all stationary solutions occur for a± = a^. 
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Figure 3. Variation of the eccentricity of each planet with Aro (left frame) 
and with cr (right frame). Initial conditions were chosen inside the QS region 
following line A in Fig.[T]for e = 0.40. The radial distance is the value of 
the osculating eccentricity. The position of the ACR solution is shown in 
red, and is located at A to = 180°. Note, however, the existence of large- 
amplitude solutions around Ato = 0, even though no stable ACR solution 
is found in this region. 

have called them Anti-Lagrangian solutions and they are connected 
to the classical equilateral Lagrangian solutions by the a-family of 
periodic orbits. By analogy, we have denoted the new solutions as: 

AL 4 : a G [0,180°] Azu g [180°, 360°] (1) 

AL 5 : a G [180° , 360°] Azu G [0°, 180°]. 

As with all previous stationary solutions, these asymmetric points 
are found at ai = a^. 

As the eccentricities grow (e.g. = 0.40) the QS region at 
(cr, Azu) = (0, 180°) causes a distortion and compression of the 
stable asymmetric domain. The Anti-Lagrangian zone seems less 
affected and surrounded by a larger island of stable motion. This 
effect is even more pronounced for = 0.60 and = 0.70 where 
the stable domain around L4 and L5 almost disappear. The region 
around AL4 and AL5 are still visible, although they also decrease 
in size and their location approaches the unstable symmetric peri- 
odic orbit located at a = Azu = 180° . 

The decrease in the size of the stable regions around the asym- 
metric ACR solutions is accompanied by a significant increase in 
the stable domains around QS orbits, which, for high eccentricities, 
seem to cover a large proportion of the plane. Inside this region we 
also note two families of periodic orbits; the Aw-family which is 
restricted to a small region around Azu = 180°, and a smaller o- 
family close to zero value of the resonant angle. 

TableQ]summarizes the detected stable stationary solutions in 
the planar planetary three-body problem, as well as their location 
in the plane of angular variables for low eccentricities. 
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Figure 4. Variation of the eccentricity of each planet (e;) with Ara (left 
frame) and with <r (right frame) for initial conditions inside the stable region 
connecting L4 and AL4 (selected from line B in Fig.[T]for e = 0.40). The 
radial distance is the value of the osculating eccentricity. The resonant angle 
Ara oscillates around one of two possible centers. One corresponds to the 
L4 configuration while the other to AL4 configuration. Locations of the 
ACR solutions are given by the fixed points in red. 



2.1 Motion Around the Stationary Solutions 

In order to visualize the dynamics of stable orbits outside the ACR, 
we integrated several orbits with initial elements 01 = 02 = 1 
AU, ei = e2 = 0.4, a = and different values of Azu. Each 
initial condition was chosen along line A drawn in Fig. Q] for e» = 
0.40. Results are shown in Figure [3] The left-hand plot shows the 
orbital evolution in the (e.2, Azu) plane, while the right-hand graph 
presents the variation of (e2, <r). In both cases the numerical output 
was filtered to eliminate short-period variations associated to the 
mean anomalies of both planets. Note that all trajectories display 
small-amplitude oscillations in a, consistent with starting positions 
near the a-family of periodic orbits. 

The behavior in the (e2, Azu) plane is more intriguing. Initial 
conditions with Azu 6 (90° , 270° ) exhibit oscillations of differ- 
ent amplitudes around the ACR solutions corresponding to quasi- 
satellite motion. Recall that this ACR solutions is located at Azu = 
180°. However, initial conditions with -90° < Azu < 90° dis- 
play regular motion that seem associated to large amplitude oscil- 
lations around Azu = 0, even though this is an unstable point lead- 
ing to close encounters and a collision between both planets. Nev- 
ertheless, there appears to be a minimum allowed amplitude for 
these solutions (shown in Figure [3] as a blue dashed curve), which 
corresponds to a semi-amplitude in Azu of approximately ~ 45° . 
Smaller amplitudes are unstable and lead to the ejection of one of 
the planets in short timescales. 

Figure [4] shows results for initial conditions inside the stable 
region connecting L4 and AL4. Semimajor axes and eccentricities 
were the same as in the previous plots. The initial values of Azu 
were varied from zero to 360° , and in each case a was chosen along 
line B in Fig.Q]for e; = 0.4 (a-family). 

The (e2, Azu) plane (left frame) shows two centers of oscil- 
lation, one corresponding to each ACR, and identified by red dots. 
L4 is located at Azu = 60° while AL4, roughly at Azu — 240° . 
As before, we see a smooth transition in the dynamical behavior 
between both modes, with no evidence of any separatrix. Conse- 
quently, it appears that any initial condition will lead to a stable 
oscillation of Azu around the nearest stationary solution. 

The motion of the resonant angle a (right frame) shows a 
different behavior. Only initial conditions very close to either L4 
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Figure 5. Orbit configuration with initial conditions chosen inside the re- 
gions of stable motions (see top left-hand frame). Initial position of planets 
are marked with open circles while crosses indicate the orbital configura- 
tion position leading to a minimum distance between the bodies. For QS, 
L4 and Z/5, the minimum distance coincides with the initial condition. 



or AL4 will show a small-amplitude circulation around the corre- 
sponding stationary point. As an example, notice some trajectories 
oscillating around a — 90° without reaching the fixed points. Fi- 
nally, due to the intrinsic symmetry in co-orbital motion, the same 
behavior is also noted for initial conditions between L5 and AL5. 

To better visualize each stable configuration, Figure[5]presents 
the orbit scheme for five stable solutions, whose initial values of the 
angles are shown in the top left-hand frame. Five initial conditions 
correspond to the stable ACR solution discussed previously (QS, 
L4, Z/5, AL4, AL5). Each of the other plots show the orbital repre- 
sentation of each solutions in (x, y) astrocentric cartesian coordi- 
nates. Initial conditions for both planets are shown in open circles, 
with mi located along the 2>axis. Both axis directions are fixed. 
The orbital trajectory of each planet (over one period) is drawn in 
thin black lines, and the configuration leading to a maximum ap- 
proach between both planets is shown with crosses. For QS, £4 
and L5, the minimum distance coincides with the initial condition. 

For QS orbits, the relative position of is always located in 
the positive x-axis, similar to the behavior noted in the restricted 
three-body problem (Mikkola et al. 2006). The relative motion of 
all five ACR solutions are periodic orbits, and symmetric with re- 
spect to the x-axis. 
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Figure 6. Families of stable QS stationary solutions in the (ei, ey) plane, 
for three different mass ratios ni2 /mi. Notice a locus of solutions at e\ = 
e2 = 0.565 for all masses. The solutions for m2/mi < 1 are mirror 
images of those for m2/m\ > 1. 



3 SEMI ANALYTICAL MODEL 

One drawback in the previous numerical approach is the excessive 
CPU time required for the construction of each dynamical map. 
In order to extended these results to other values of the parameter 
space (e.g. planetary masses, eccentricities) it is useful to construct 
a semi-analytical model for the co-orbital motion. 

Such a model can be developed along similar lines as other 
mean-motion resonances (e.g. Michtchenko et al, 2006, 2008a, b). 
It requires two main steps: first, a transformation to adequate res- 
onant variables and, second, a numerical averaging of the Hamil- 
tonian with respect to short-period terms. Both tasks are detailed 
below. 

We begin introducing the usual mass-weighted Poincare 
canonical variables (e.g Laskar 1990) for each planet m,: 



A, 
—, 




(2) 



where fii — c 2 (mo + mi), k denotes the gravitational constant, 
and m'i is the reduced mass of each body, given by: 

/ m 4 m 

TO, = ; . (3) 

n%i + mo 

The Hamiltonian function F can be expressed as F — Fa + Fi, 
where Fo corresponds to the two-body contribution, and has the 
form: 



Fo 



2L? ' 



(4) 



For initial conditions in the vicinity of co-orbital motion, we 
define the following set of planar resonant canonical variables 
(Ji, I2, IC, AM, a, Avj, Q, q), where: 



o — A2 — Ai 

Azu = VJ2 — V3\ 
q = ZU2 + ZUl 

Q — Ai + A2 — q 



h 
h 
Jl 

J2 = |(L 



1(^2- 

i(G 2 - 



Li) 

Gi- 
-G 2 ) 
■L2) 



L2 + L1 



(6) 



where, Ji = ^AM and J 2 = i/C. A generic argument ip of the 
disturbing function can be written as: 



<P = jl Al + J2A2 + js^l + jiZU 2 



(7) 



where jk are integers. In terms of the new angles the same argument 
may be written as: 



Kh ~ h)<r + Ui - h)Azu + (ji + j 2 )Q] ■ 



(8) 



Since q is a cyclic angle, the associated action AA4 is a constant of 
motion (total angular momentum) of the system. 

The next step is an averaging of the Hamiltonian over the fast 
angle Q. This procedure can be performed numerically, allowing to 
evaluate the averaged Hamiltonian F as: 



F{h,h,a, Aw, K, AM) 



1 

2^ 



FdQ. 



(9) 



In the averaged variables, IC is a new integral of motion which, in 
analogy to other mean-motion resonances (e.g. Michtchenko et al. 
2008a), we call the scaling parameter. 

F then constitutes a system with two degrees of freedom in the 
canonical variables l2,o, Azu), parametrized by the values of 
both AC and AM. Since the numerical integration depicted in equa- 
tion (|9j is equivalent to a first-order averaging of the Hamiltonian 
function (e.g. Ferraz-Mello, S. 2007), only those periodic terms ifTJ 
with ji + j% = remain in F. In consequence, we can rewrite the 
generic resonant argument of the averaged system as: 



(10) 



where the index 32 , ji are integers that may take any value in the 
interval (— 00, 00). 



4 FAMILIES OF PERIODIC ORBITS 

In the averaged system defined by F exact zero- amplitude ACR 
solutions are given by the stationary conditions: 



dF 
9^ 



dF 
dAvo 



dF 
dh 



dF 
W2 



0, 



(11) 



and can therefore be identified as extrema of the averaged Hamil- 
tonian function. In the present section we will use this approach to 
estimate the families of different ACR as function of the planetary 
masses and eccentricities, and compare the results with numerical 
integrations of the exact equations of motion. 



The second term, F± , is the disturbing function which can be writ- 
ten as: 



2 1 
Fi = —ft mim 2 — + Ti, 



(5) 



where A is the instantaneous distance between the two planets, and 
Ti is the indirect part of the potential energy of the gravitational 
interaction (see Laskar 1990, Laskar and Robutel 1995 for more 
details). 



4.1 Families of Symmetric ACR. QS 

We begin calculating the exact stationary solutions, correspond- 
ing to QS configurations, as a function of the eccentricities, and 
for different values of the planetary masses. As mentioned in Had- 
jidemetriou et al. (2009), the locations and stability of the ACR do 
not appear dependent on the individual values of the masses, but 
only on their ratio 7712 /mi . 
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Figure 7. Equilibrium values of a and Am for the family of AL4 solutions 
as function of the eccentricity of the smaller planet, for several mass ratios 
mz/mi > 1. 




Figure 8. Families of AL4 solutions, e2 as function of ei, for several mass 
ratios 1712/ mi > 1. 



In all cases, the stationary values of the canonical momenta 
Li are such that n\ = 712, where rii are the mean motions of 
the planets. For equal mass planets, this reduces to the condition 
ai = a,2. Finally, the angles of the exact ACR always remain 
locked at (a, Aw) = (0, 180°). Hadjidemetriou et al (2009) pre- 
sented similar plots for the same mass ratios. 

Figure [6] shows the families of stable zero-amplitude QS or- 
bits for selected mass ratios: 7712/mi = 1/3, 7712/7711 = 1, and 
milm\ > 1. For equal masses, all solutions occur for and e 1 — e^. 
Due to the intrinsic symmetry of the dynamical system, the family 



of stationary solutions for 7712/ 'mi = 1/3 is a mirror image of the 
solution for m%lm\ = 3, since it may be obtained by simply in- 
terchanging ei with e2. In the case of 7712/mi — 3, we note that 
e2 < ei for e2 < 0.565, while e-2 > ei for more elliptic orbits. 
Figure [6] shows also the solutions for 7712/777.1 = 2, 5, 20 mass 
ratios. For mass ratios smaller than unity, the solutions are mirror 
images with respect to the family 7712 /mi = 1. Note that the fami- 
lies of stable solutions approach ei = e2 as 7712 — » mi . However, 
as the mass ratio tends towards the restricted three-body problem, 
the eccentricity of the smaller mass approaches unity. Finally, the 
solution a — 0.565 is common to all the QS families, and cor- 
responds to a global extrema of the Hamiltonian in this plane. A 
similar structure was already noted by Michtchenko et al. (2006) 
for other mean-motion resonances. 



4.2 Families of Asymmetric ACR Solutions. L4 and AL4 

The same procedure can also be applied to the Lagrangian L4 and 
Anti-Lagrangian AL4 configurations. Recall that the dynamical 
maps (Figure [2} showed a symmetry with respect to the transfor- 
mation (a, Azu) — > (— a, —Aw), so the results discussed here can 
also be applied to the L5 and AL$ solution, by applying the same 
operation on the variables. 

The ACR solution associated to the Lagrangian solution L4 
shows no variation in the angles, maintaining constant both angles 
at 60°. The solutions remain stable for initial conditions up to ec- 
centricities ei — 0.7. However, the AL4 shows significant changes 
as function of the eccentricities. Figure [7] shows the equilibrium 
values of both angles for the family of AL4, as a function of the 
eccentricity of the smallest planet, for several values of the mass 
ratio 77i2 /mi. The resonant angle a increases monotonically from 
60°, at quasi-circular orbits, towards ~ 180° for near parabolic tra- 
jectories. As the mass ratios increases, the maximum value of the 
resonant angle decreases, reaching a = 150° for a mass ratio of 
m%lm\ = 10. 

The secular angle Azu shows a slightly more complex behav- 
ior. Initially it increases from ~ 240° until it reaches a maximum 
value close to ~ 260°, after which it once again decreases towards 
Am ~ 180°. The planetary eccentricity corresponding to the max- 
imum in the secular angle increases with the mass ratio, approach- 
ing the parabolic limit for m^i 'mi ~ 10. 

As shown in Figures [Tj and [2] the size of the stable region 
around each asymmetric solution decreases with the increase of d, 
and practically disappears as the angles approach 180° degrees. For 
quasi-parabolic orbits, only the region around AL4 is discernible. 
Thus, for high eccentricity planets in co-orbital motion, it appears 
that the AL4 and AL$ asymmetric solutions are more regular than 
the classical equilibrium Lagrangian solutions L4 and L5. 

The values of the planetary eccentricities at AL4 for different 
mass ratios is presented in Figure[8] Contrary to the QS trajectories, 
there appears to be a purely linear dependence between 62 and e\ 
as a function of the mass ratio. In fact, a simple numerical analysis 
of the results appears to indicate that 



ei 



7712 

mi 



("2- 



(12) 



Thus, for mass ratios approaching the restricted three-body prob- 
lem (with 7712 — > 0) it should be expected that the eccentricity 
of the massive planet mi at the AL4 solution would tend towards 
zero. 

Finally, the equilibrium values of the semimajor axes also 
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Figure 9. Exact numerical integrations of initial conditions close to the 
AL4 stationary solutions (m.2/ 'mi > 1). Black curves represent the AL4- 
family of ACR calculated with the semi-analytical model. 



change as function of the mass ratio. Here, however, it is easy to 
see from the stationary conditions i ll It that a zero- amplitude AL4 
trajectory is characterized by the relation m = 712. For equal mass 
planets, this reduces to ai = 02- 

The families of stationary solutions presented in this section 
were calculated using our semi-analytical model. In order to com- 
pare them with actual numerical simulations of the exact equa- 
tions, we choose four solutions from Figure [7] with ei = 0.2, 
but corresponding to different mass ratios. Each was then numer- 
ically integrated for several orbital periods, assuming zero ini- 
tial values for the cyclic angular variables q and Q. Results are 
shown in Figure [9] where the top frame presents the trajectories 
in the plane (ei cos a, ei sin a) and the bottom graph in the plane 
(e2 cos A137, e2 sin Avj). Each initial condition shows a small am- 
plitude oscillation around the stationary value, which presents a 
very good agreement with the family of AL4 solutions calculated 
with our model (black curve). 



4.3 Adiabatic Mass Variation in AL4 

As a final analysis, in this section we study the orbital evolution of 
a system initially near AL4, when the mass of one of the planets 
is decreased adiabatically. This question is raised for three reasons. 
First, as shown by Lee (2004), for two planets in a 2/1 mean-motion 
resonance, a sufficiently slow change in one of the masses will pre- 
serve the resonant configuration and allow to calculate the variation 
of the ACR as a function of mi/mi. In other words, this approach 
provides a different numerical test of our semi-analytical model and 
an alternative way to calculate the stationary orbits. Second, the re- 
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Figure 10. Evolution of AL4 applying a smooth adiabatic decrease in mi. 
Initial conditions correspond to m,2 = mi and e\ = ei = 0.2. The 
stationary values of the angles are invariant to the mass change, although 
the amplitude of Aro increases as e2 approaches zero. The ratio of the 
equilibrium eccentricities follow the relation U2t . as shown by the dashed 
black curve overlaying the data of C2- 



suits will also allow us to test the robustness of the new asymmetric 
co-orbital solutions AL4 and see how they respond to changes in 
the parameters of the system. Finally, we wish to analyze the behav- 
ior of these new solutions in the limit of the restricted three-body 
problem, corresponding to mi = 0. 

Figure [JjJJ shows a typical example. Initial conditions corre- 
spond to an AL4 solution for m.2/mi = 1 and ei = 0.2. While 
mi was maintained fixed, was varied linearly down to mi = 
in a timescale of 10 6 orbital periods. We checked using other 
timescales, finding no significant variations. This guarantees that 
we are effectively in the adiabatic regime. 

The top graph of Figure [lOl shows the evolution of the orbital 
eccentricities as function of the mass ratio. As soon as 7712/mi 
departs from unity, the value of e2 increases while e\ decreases. 
The broken black curve that can be seen over the red curve shows 
the predicted value of ej applying the relation d 1 ~2 b to each value of 
62. The agreement is excellent, giving an additional corroboration 
to this empirical relationship between the eccentricities. It must be 
noted that neither the total angular momentum AM nor the scaling 
parameter K, are preserved during the mass change. The bottom plot 
of Figure [TO] shows the behavior of the angular values during the 
mass variation. The equilibrium values of both a and Azu remain 
practically unchanged. 

For smaller mass ratios mi tends towards a circular orbit, 
while the eccentricity of the smaller planet approaches e2 ~ 0.4. 
This seems to imply that the asymmetric AL4 (and consequently 
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Figure 11. Orbital trajectories in a cartesian rotating pulsating reference 
frame, where the positions of mo and mi are fixed in the x-axis with unit 
mutual distance. Both bodies are marked with large filled circles. The brown 
circle shows the orbital evolution of ni2 when placed in L4, while red dots 
correspond to an initial condition for m,2 placed in AL4 (see Figure fTol . In 
both cases the plot presents the orbital evolution as m,2 — > 0. No change 
is observed in L4, while the blue curve shows the final orbital trajectory 
around AL4 when ni2 reaches zero. The oscillation period is equal to the 
orbital period between the massive primaries. 



AL$) solutions could also exist in the limit of the restricted three- 
body problem. To test this conjecture and compare the trajectories 
of both Z/4 and AL4 solutions in the restricted (m.2 — > 0) limit, 
Figure [Til plots the (x, y) cartesian coordinates of two initial con- 
ditions in a rotating pulsating reference frame. 

In the rotating pulsating reference frame the positions of both 
mo and mi are fixed in the a;-axis (shown with large solid circles). 
Three orbital evolutions are shown: the brown dot corresponds to 
initial conditions in the asymmetric L4, while red dots map the evo- 
lution of an orbit originally in AL4. In both cases we started with 
7?i2 = mi, but subsequently decreased m,2 to zero (restricted case). 
No change is noticed in the L4 orbit, and the trajectory remained 
in an equilateral configuration with the two finite masses. However, 
the AL4 solution converged towards a tadpole-type orbit of large 
amplitude (blue curve) for 7712 — > 0. This solution corresponds to a 
periodic orbit whose period coincides with the orbital period of the 
primaries around the center of mass. Green dots map the evolution 
of an orbit originally in QS. As we can see the orbit described by 
QS configuration revolves around the mi, in the same way that was 
observed in the restricted problem. 

Thus, there appears to be a structural difference between the 
Z/4 and AL4 planetary solutions discussed in this paper. Although 
both appear as ACR (fixed points in the averaged problem) the first 
are true stationary solutions in the unaveraged rotating frame, while 
the new solutions AL4 are actually large amplitude periodic orbits 
that encompass the classical Lagrangian equilateral solution. 



5 CONCLUSIONS 

We studied the stability regions and families of periodic orbits 
of two-planet systems in the vicinity of a 1/1 mean-motion reso- 
nance (i.e. co-orbital configuration). We considered different ratios 
of planetary masses and orbital eccentricities, also we assumed that 
both planets share the same orbital plane (coplanar motion). 

As result we identified two separate regions of stability, each 
with two distinct modes of motion: 

• Quasi-Satellite region: Originally identified by Had- 
jidemetriou et al. (2009) for the planetary problem, QS orbits cor- 
respond to oscillations around an ACR located at (a, Avj) = 
(0, 180°). Although not present for quasi-circular trajectories, they 
fill a considerable portion of the phase space in the case of moderate 
to high eccentricities. 

We also found a new regime, associated to stable orbits display- 
ing oscillations around (a, Azu) = (0, 0), even though this point is 
unstable and corresponds to a collision between the two planets. 

• Lagrangian region: Apart from the previous symmetric solu- 
tions, we also found two distinct types of asymmetric ACR orbits 
in which both a and Aw oscillate around values different from 
or 180°. The first is the classical equilateral Lagrangian solution 
associated to local maxima of the averaged Hamiltonian function. 
Independently of the mass ratio mj/mi and their eccentricities, 
these solutions are always located at {a, Aw) = (±60°, ±60°). 
However, the size of the stable domain decreases rapidly for in- 
creasing eccentricities, being practically undetectable for > 0.7. 

The second type of asymmetric ACR correspond to local min- 
ima of the averaged Hamiltonian function. We have dubbed them 
Anti-Lagrangian solutions (AL4 and AL5). For low eccentricities, 
they are located at (a, Am) — (±60°, =Fl20°). Each is connected 
to the classical L4 and L$ solution through the cr-family of periodic 
orbits in the averaged system. Contrary to the classical equilateral 
Lagrangian solution, their location in the plane (a, Azu) varies with 
the planetary mass ratio and eccentricities. Although their stabil- 
ity domain also shrinks for increasing values of ej they do so at 
a slower rate than the classical Lagrangian solutions, and are still 
appreciable for eccentricities as high as ~ 0.7. 

Finally, we also applied an ad-hoc adiabatically slow mass 
variation to one of the planetary bodies, and analyzed its effect 
on the AL4 configuration. We found that the resonant co-orbital 
solution was preserved, with practically no change in the equilib- 
rium values of the angles. The eccentricities, however, varied with 
the larger planet approaching a quasi-circular orbit as the smaller 
planet had its eccentricity increased. These solution still exist in the 
limit of the restricted three-body problem (i.e. mi — > 0), although 
both types of asymmetric solutions (L4 and AL4) have different 
geometries. While the first are true stationary solutions in the un- 
averaged system, the latter are periodic orbits around the classical 
equilateral Lagrangian points. 
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